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Abstract 

We study the critical behavior of the conductivity a{u) at the zero tempera- 
ture superconductor-Mott insulator transition in d space-time dimensions for 
a model of bosons with short-range interaction and no disorder. We obtain 
a{<jJn) = (4e 2 /ft)<7 (E a4 -<E , as predicted by the scaling theory, and the pref- 
actor a e is calculated in the e-expansion, to order e 2 (e = 4 — d). In two 
spatial dimensions, (d = 3), we find a value of the universal conductance 
a = 0.315(4e 2 //i), in good agreement with the known Monte Carlo results. 

PACS numbers: 74.75 +t, 74.40 +k 
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In two dimensions at the T = superconductor-insulator (SI) transition the conduc- 
tance is finite and universal^!. Since the original prediction of a metallic behavior at zero 
temperature for two dimensional superconductors there has been a lot of interest both on 
experimental and theoretical side. 

The systems where the T = SI transition has been observed are ultra thin super- 
conducting films! and Josephson junction arrays@~i. The films are driven through the SI 
transition by varying the thickness while, for the arrays, the ratio of the charging energy to 
the Josephson coupling determines whether they are in the insulating or in the supercon- 
ducting phase. 

Many theoretical works have investigated various aspects of the SI transition. The 
superconductor-Mott insulator transition and the universal conductance in a model with 
no disorder was considered in Ref. by means of 1/N expansion and Monte Carlo simu- 
lations. The duality between charges and vortices, leading to an universal conductance of 
(4e 2 //i)@, has been discussed in Ref. [| Fisher et a/.i studied the disordered interacting boson 
problem using a scaling approach. The dirty boson system and the transition to the Bose 
glassphase (including the case of long-range Coulomb interaction) was extensively studied 
in0H Wen employed a scaling theory of conserved current at anisotropic critical points@ 
identifying many universal amplitudes. One of these amplitudes in two dimensions reduces 
to the universal conductance a* . The finite frequencies properties close to the transition 
point were analyzed by means of the 1/N expansionE10. 

The determination of the value of a* relies almost entirely on numerical methods (Monte 
Carloi and exact diagonalizationlll). The only analytical calculation of a* is based on 
the 1/N expansion! Another powerful and widely used method to calculate the critical 



properties is the e-expansion (for a review see Ref. 17). 

In this Letter the conductivity at the SI transition is calculated applying the e-expansion. 
We will compute the scaling form of the cr{uS) discussed in Refs. |H,|l2l and we will determine, 
in two spatial dimensions, the value of a*. 

The description of the SI transition starts from the hypothesis that the critical properties 
are determined by the quantum fluctuations of the phase of the superconducting order 
parameter. A widely used model for the SI transition is a Hubbard like model for strongly 
interacting bosons on a lattice (Bose-Hubbard model) 

i i (ij) 

where b\b are the creation and annihilation operators for bosons and rij is the number of 
bosons, [/describes the short range interaction between bosons, \x is the chemical potential 
and t is the hopping matrix element. At integer fillings, the Bose-Hubbard model in d — 1 
spatial dimensions is in the same universality class of the XY model in d dimensions. It 
is more convenient for computational purposes to work with the Ginzburg-Landau action 
(onto which the Bose-Hubbard model maps) 

F[1>\ = J d^rdr {<W| 2 + | W| 2 + \d T ij\ 2 + ^H 4 } • (2) 

where ip(r, r) is a complex one-component scalar field. In (0) we have introduced e = 4 — d 
and the arbitrary momentum scale k in order to keep the quartic coupling m dimensionless. 
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The bare parameters So and uq can be related to the coupling constants in the Bose-Hubbard 
modeli'El. Mo includes the effect of the interaction and 5q measures of the distance from the 
critical point (we set Sq = in the following). It is evident from eq. (|2|) that the quantum 
system at zero temperature is equivalent to a classical system in d dimensions. 

The conductivity of the system can be obtained by minimally coupling the scalar field 
in eq. (0) to a vector potential A. The imaginary time frequency-dependent conductivity is 
given by 



h 



d d x rdr 



S 2 InZ 



5A„(r,r)5A v (0) 



(3) 



A=0 



where Z is the partition function of the system. The conductivity is obtained by analytically 
continuing eq. @ to real frequencies (icu n = u + i0 + ). Being interested in the zero field 
case, we will only consider the longitudinal conductance (e.g. a xx {oj n ) ), which we shall 
generically indicate as a(u n )). 

By performing the functional derivatives it is straightforward to express the conductivity 
in terms of two and four point Green's functions. Following Cha et al.Q it is convenient to 
introduce the notation a{u n ) = (4e 2 /h)p(u n )/uj n where 



d d q 
(2tt)° 



d d q d d p 
(2vr) d (2n) d 



p+ 



(4) 



In eq. (||) the q are real vectors in the d- dimensional space and k = (0, 0, ...,u n ). 

The two-point function in eq. (§), expressed in terms of the self-energy T,(q) is (in the 
following we shall use the notation q = \q\): (V$0<f) = G{q) = (q 2 + S(g)) x . Since the three 
point irreducible vertex is zero in the insulating phase where no spontaneous symmetry 
breaking occurs, the four point function in eq. (f|) is the sum of a connected part containing 
the irreducible four point vertex Y^\q,p, k), and two non-connected parts of which only one 
gives non vanishing contribution to p(uj n ). T^(q,p, 0) and £(g) are related by the following 
Ward identity 



dq x 



(2tt 
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This identity, related to the underlying gauge symmetry of the free energy in eq.(^), is 
essential to prove that p{oJ n = 0) = 0. This can be easily checked by evaluating the derivative 
of the two-point Green's function and substituting in eq. (ffl). In the superfluid phase the 
cancelation does not occur due to the presence of three point vertex in the theory. 

In order to get a more suitable expression of p{u n ) for a perturbative analysis, we split 
the four point irreducible vertex in five parts, T^(q,p, k) = + F^(q + p) + f\ (q — p) + 
4 r)+^res(9)Pi ^) according to the external momenta dependence. has no momentum 
dependence and T^ s (q,p, k) is the residual part with momentum dependence different from 
q + p or q — p or k. This splitting turns out to be interesting because and r^ 4 )(/c), when 
inserted in eq. ([|), vanish and F$ s (q,p, k) gives a contribution to eq. (f|), in the e-expansion, 
only of order 0(e 3 ). Eq.(0) and the decomposition of the I^ 4 ) allow, after some algebra, to 
express p{uJ n ) in the form 
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d d q d d p 



oq x 



G{q) G{q)-G{q + k) 



J ^-^ ) - d q x p x G(q)G(p})^[G(q + k) - G(q}\ 7 W (<? P) [G(p+ k) - G(p) 
-G(fi [2 G(q + k) - G(q)} rH(q,p,0) + G(p+ k)G(q + k)4l(q + |, k) 



(6) 



where we have introduced "f^(q — p) = 1$ (q '— p) — T^(q — p). So far no approximations 
have been used to obtain eq. ([]). The next step is the computation of the self-energy and 
four point vertex in the framework of the e -expansion. 

At the critical point G~ l (q) = q 2 ~' n and the value of the critical exponent rj, for the 
theory considered, is known to be 77 = (e 2 /50)(l + 19e/20) + 0(e 4 ). The expansion of G(q) 
in powers of e follows straightforwardly. 

Let us now consider the four point vertex T^'(q,p, k); its expansion in terms of the 
dimensionless renormalized coupling constant u can be formally written as 



r (4) («f,p, k) = K e \-2u + u 2 A(q,p, k) + u 3 B(q,p, k) + 0(u 4 



(7) 



It is understood that the function in eq. (0) is renormalized and therefore, employing 
dimensional regularization, no pole at e = is present in A(q,p,k) and B(q,p,k). In the 
following calculation the minimal subtraction renormalization prescription has been used. 

Since we are concerned with the computation of the four point vertex at the transition, 
we must evaluate eq. (0) at u = u*, u* being the fixed point of the (3 function of the 
theory: (3{u*) = 0. By introducing the counterterms ai through the expansion of the 
bare coupling w appearing in eq. (0): u = u(l + a\U + a 2 u 2 + ...), we find from the 
evaluation of the bare four point vertex at two loop in the minimal subtraction scheme: 
ai = 5/(8vr 2 e), a 2 = (25/e 2 - 15/(2e)) /(8vr 2 ) 2 . 

The fixed point u*, obtained as a function of the counterterms dj (see for instance Ref. 
Mj, is u* = (8tt 2 ) (e/5 + 3e 2 /25) + 0(e 3 ) . Since the coefficients A(q,p, k) and B(q,p, k) do 



not contain any pole at e = 0, a full 0(e 2 ) calculation of the four point vertex requires the 
computation of the one loop part A(q,p, k) only. 

Nevertheless, let us consider the momentum dependent structure of both A(q,p,k) 
and B(q,p,k). In the one loop part A(q,p,k) one has terms like e n log n+1 , e n+1 log n+1 , 



r-n+2 



log 



n+l 



with n = 0, 1,2, .. and log being the logarithm of a generic combination of 



external momenta; in the two loop part B(q, p, k) one gets e n log n+2 , e n 



log 



n+l ,n+l 



log 



n+l 



As mentioned above, a full 0(e 2 ) computation of eq. (|?p requires just the coefficient of the 
term e° log in A(q,p,k), but, since we are interested in summing formally the leading log 
series in eq. (|7|), we need to know at least the first two terms of the series. This amounts to 
calculate the terms e log 2 in A(q,p, k) and e° log 2 in B(q,p, k). 
In the minimal subtraction scheme A(q,p, k) is 



— * 

A(q,p,k) 



c 



1 W+p\ , \q-P\ 

log 2 log 



K 



\k\ 

2 log^ 

K 



+- log 2 ^ — ^ + e log 2 ^ — ^ + e log 2 i-L + 0(e) 

L, Kj Kj Kj 



(8) 



4 



where, according to our task, we have displayed all the contributions which are relevant for 
the computation of the first two terms of the leading log series and neglected the subleading 
logarithms (e log) and the higher powers of e. C in eq. @ is a 0(1) number, irrelevant for 
the determination of p(uj n ) in eq. (H). 

Still we need the e° log 2 part of B(q,p, k), which we indicate as B(q,p, k); it is 



B(q,p,k) 
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2\2 



log 2 ^ + 2 log 2 ^ 

K 



(9) 



As stated above, to order e 2 , r^ s (g, p, k) = 0, making eq. (|6|) much simpler. Actually 
the first nonvanishing contribution to T^ s {q, p, k) comes from the subleading term e° log 
in B(q,p,k) which is neglected in our calculation. We are now able to sum formally the 
leading log series obtaining some power of the momenta involved (note that the momentum- 
independent term 0(e) in the series is irrelevant because it vanishes when integrated in eq. 



(0)). Putting together eqs.([7|, 



7^ (<l~p) becomes 



7 (4) (? -P) 



167r 2 e 
25 



g- P T + 0(e 3 ) 



(10) 



The importance of summing the leading log series is now clear: eq. ( |T0D is correct only to 
order 0(e 2 ), but the infrared singularities due to the logarithms in eqs. (||, |9|), disappear in 
the resummed quantity \q — pf. 

Another comment is in order: the last result does not depend on the arbitrary momentum 
scale k introduced in eq. nor on the choice of the minimal subtraction scheme. Indeed, 
one can easily prove that the shift of the counterterms by a finite (no poles at e = 0) 
amount, <2j — > Oj + Hi, does not change 7 ( - 4 - ) ((f — p) in eq. (|10|). 

Finally, going back to eq. (|6|), we can insert eq. (|TUD, and discard the terms containing 
r^ s (g, p, k), thus obtaining p(cu n ) and, consequently, a(cu n ) to order e 2 ; in dimensions lower 
than 4 all integrals in eq. (|^) are convergent in the ultraviolet (in d > 4, a short-range cutoff 
should be introduced to regularize the results, but this is not important for our purposes) 
and we find 




where 

hi 
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dX 



r 2 (d/2)r(i - d/2) 
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(— a(l — a)\ — iay^j 



and z(d) 



-i 



(11) 



(12) 



T(x) and 



\ d ' 2 /{T{d/2)) J °° dx x d ' 5 (l - {x/2f- d/2 T (d/2) J m -i)(x) 
J v (x) are, respectively, the Gamma and Bessel functions and the u n independent factor 
in eq. ( |TTD has been expressed in integral form for the sake of simplicity. 

Eq. ([TT|) is the central result of our paper. The form of the frequency dependent conduc- 
tivity, as predicted from the scaling theoryEE, is obtained. 

In the calculation of the a(u n ), e-expansion was used to determine the irreducible vertices. 
The d- dimensional integrals in eq.(^) have been instead performed not resorting to any 
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expansion. This observation is crucial in order not to obtain spurious dependences on the 
ultraviolet cut-off needed to regularize the theory in higher dimensions. In this sense the 
calculation of the conductivity to 0(e 2 ) should be understood as the scheme of approximation 
for the vertices. 

In two spatial dimension the universal conductance a* can be evaluated from eq. ( |TT|) 
at d = 3; this leads to the result displayed in Table 1, together with the other known 
determinations of cr*. Our estimate is in agreement with the results of the Monte Carlo 
simulations in Ref. [7] within 10% but it should be noted that the 1/N and e-expansion 
approach the Monte Carlo value from different sides. 

As a final remark we want to go back to the connection with the experiments. The 
model with only phase fluctuations is certainly applicable to Josephson junction arrays, but 
it has been recently questioned in the case of thin films. New measurements^ point out 
that the value of the conductance may be not universal and this should be associated to 
the relevance of fermionic degrees of freedom at the transition. This is also the outcome of 
recent calculations which include fermions starting from a Bose model with local damping0 
or simulating a model of electrons interacting through an attractive potential^. 
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TABLES 



Method 




(l/iV)-exp w 




= 


0.315 


0.251 


0.285 



TABLE I. The value of the universal conductance obtained by various methods are reported 
for comparison, a) this work, b) - c) Ref. ^ 
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